155 1561-5359. Штучний інтелект, 2017, Хо 1 


УДК 536.421 


М.Г. Бердник 
Державний вищий навчальний заклад «Національний гірничий університет», Україна 
пр. Карла Маркса, 19, м. Дніпро, 49005 


МАТЕМАТИЧНА МОДЕЛЬ І МЕТОД РІШЕННЯ 
УЗАГАЛЬНЕНОЇ ЗАДАЧІ ДІРІХЛЕ ТЕПЛООБМІНУ ТІЛА, 
ЯКЕ ОБЕРТАЄТЬСЯ 


М.С. Вегапук 
Утаїе Нієрег Едисайоп Пяціицоп "Мацпопа! Міпіпеє Юпімег8ісу", ОКгаїпе 
19, Кагі! Магх ау., Опірго, 49005 


МАТНЕМАТІСАЇ, МОРЕІ, АХЬр СЕХЕКАТЛЛЕР 5ОГОТІОМ 
МЕТНОЮР РІКСНЕЇ.Е РКОВІ.ЕМ5 ОЕ НЕАТ ЕХСНАХСЕ ОЕ 
КАРІАТТОХ ВОрУ 


У статті за допомогою розробленого нового інтегрального перетворення знайдено 
температурне поле ізотропного тіла обертання з відомим рівнянням твірної лінії, яке обмежене двома 
торцями і бічною поверхнею тіла, яке обертається з постійною кутовою швидкістю навколо осі 07, з 
урахуванням кінцевої швидкості поширення тепла, у вигляді збіжних ортогональних рядів по 
функціях Фур'є. 

Ключові слова: тіло обертання, інтегральне перетворення, критерій Фур'є, час релаксації. 

І Фе агіїсіев, мій Ше Пер ої Ше деуеїореай пем/ іпіета! (гап5їогтацйоп, Фе (еппрегаште Пе!д ої ап 
і50їгоріс Боду ої госайоп ув а Кпоууп едцайоп ої Ше ргодисйїоп Ппе, місії 15 Пйтігга Бу (м/о епа5 апа 
Їаїега! зигвасе ої гокайоп, у/ПпісП гоіатез аї а соп5капі апеціаг уеіосісу агоцпа Ше О7 ахіз, 15 дегіуед, ур ре 
ассошиі ої їре па! гаїе ої реаї ргораваноп , іп Ше їогт ої сопуегеепі огірогопа! 5егіе5 їп Коигіег Гипсйопе. 

Кеу уогаз: Боду ої гобайоп, іпіеота! ігап5огтайоп, Боштег сгіїегіоп, геіахайоп йте. 


Вступ 

Проблема дослідження температурних полів у тілах, що обертаються, постійно 
привертає увагу дослідників, так як багато елементів машин і механізмів (супутники, 
сортопрокатні валки, ротори енергетичних агрегатів, дискові гальма та ін.) мають їх 
форму і працюють в умовах інтенсивного нагріву. Більшість робіт в теорії 
теплопровідності присвячено вивченню та аналізу температурного поля в нерухомих 
тілах обертання. У деяких роботах вивчається температурне поле в тілах при 
рухомих джерелах тепла, в інших роботах досліджується температурне поле в тілах з 
рухомими межами. 

Постановка проблеми 

При високих інтенсивних нестаціонарних процесах, що спостерігаються, 
наприклад, при вибухах, надзвукових потоках, великих швидкостях обертання вплив 
скінченності величини швидкості поширення тепла на теплообмін стає помітним П- 
2). Ось чому до числа проблем, що являють великий теоретичний і практичний 
інтерес, належить проблема вивчення температурного поля в тілах, які обертаються 
навколо своєї осі, з урахуванням скінченності величини швидкості поширення тепла. 

Аналіз останніх досліджень і публікацій 

Як показує огляд літератури, теплообмін у тілах, які обертаються, вивчений в 
даний час ще недостатньо |1|. Показано, що чисельні методи дослідження 
нестаціонарних неосесиметричних задач теплообміну циліндрів, які обертаються, є 
не завжди ефективними, якщо мова йде про обчислення при великих швидкостях 
обертання. Так доводиться |1|, що умови стійкості обчислень у методі кінцевих 
елементів і методі кінцевих різниць, що застосовуються до розрахунку 
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нестаціонарних неосесиметричних температурних полів циліндрів, які обертаються, 
визначаються аналогічними характеристиками. Ці умови мають вигляд: 


ЛЕ 
Іона вії 1 ЕН ЗР 


ЛФ! ЛФ 2 
де Бо - критерій Фур'є, Ра - критерій Предводітелева. 

Якщо Ра -10), що відповідає кутовій швидкості обертання металевого 
циліндра оаібТісек радіусом 100 мм, змінні ДФ и ДЕ, повинні бути 
підпорядковані таким умовам: 

Др «2410? і ДЕ, «2-1079, 


Для рівномірно охолоджуваного циліндра за умови Ві - 5 (Ві - критерій Біо) 
час, необхідний для того, щоб температура досягла 90905 стаціонарного стану, 


дорівнює Ко 0.025. Це означає, що потрібно принаймні здійснити 1.3.108 
операцій по часу для того, щоб було досягнуто стаціонарний розподіл температури. 
Більше того, потрібно відзначити, що протягом одного циклу обчислень 


потрібно здійснити 314.10? обчислень, так як внутрішній стан у кільці 


характеризується 3.14: 10? точками. У результаті видно, що число обчислень, 
необхідних для отримання чисельного результату, видається нереальним. Тому для 
вирішення крайової задачі, яка виникає при математичному моделюванні 
нестаціонарних процесів теплообміну в тілах, які обертаються, будемо застосовувати 
інтегральні перетворення. 

Мета дослідження 

Розробка нової узагальненої математичної моделі температурних розподілів у 
тілах, які обертаються навколо своєї осі, у вигляді крайової задачі Діріхле математичної 
фізики для рівняння теплопровідності, та розв'язання отриманої крайової задачі, 
розв'язки якої використовуються під час керування температурними полями. 

Виклад основного матеріалу 

Розглянемо розрахунок температурного поля тіла обертання (рис.1) з твірною 
лінією г - 2(2) у циліндричній системі координат (р, Ф, гай Ізотропне тіло обертання 
обмежене двома торцями 5; (2-0 ), 5» ( с -П) 1 бічною поверхнею обертання 53, 


яка перетинається з поверхнями 5; уздовж ліній /., ізі2 


Рис. 1. Тіло обертання з твірною лінією г (2). 
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Тіло обертається навколо осі О7 з постійною кутовою швидкістю и, а 
швидкість поширення тепла є відомою величиною. Теплофізичні властивості тіла не 
залежать від температури, а внутрішні джерела тепла відсутні. У початковий момент 
часу температура циліндра постійна Су, а на бічній поверхні циліндра температура 
відома і не залежить від часу У(о, 1) На торцях відомі значення температури 
СІ (г, ф) і С» (г, ф) при 7-0 і 7-В відповідно. 

В П| отримано узагальнене рівняння переносу енергії для рушійного елемента 
суцільного середовища, з урахуванням скінченності величини швидкості поширення 
тепла. Згідно з |1| , узагальнене рівняння балансу енергії твердого тіла, яке 
обертається, з постійною кутовою швидкістю (» навколо осі О7/, теплофізичні 
властивості якого не залежать від температури, а внутрішні джерела тепла відсутні, 
в циліндричній системи координат набуває вигляду: 


сто сто Годто Єт дт тост 10т даті 

ус рфФфооят, зФ вед З З Я (1), 
ді дФ ді? ОФОЇ дн? г г ? до да? | 

де у - щільність середовища; с-питома теплоємність; Т (р, Ф, 2,1) - 


температура середовища; Л - коефіцієнт теплопровідності; (- час; т,- час релаксації. 

Математично задача визначення температурного поля циліндра складається в 
інтегруванні диференціального - рівняння теплопровідності (1) в області 
,р- р, Ф, 2, др є (0, (2), Фе (0,27), 7є (0, І), (є (0, го) |і що, З урахуванням 


прийнятих допущень, запишеться у вигляді: 


90 00 | 070 9280  а|д20 100 1 070 | д2в 
чан а ОО пере важгарчюс зір Звог (2) 
ді де ді 9901  К"|др" Ро р'доф дх 
з початковими умовами 
00 ЧУ 0 
(р, р, 20)- 0, Збірирьаб) -0 (3) 
іграничними умовами 
ДСО фу сг) є (ру г), 4 
д(р, ф.0,г) - Ор, р) др.а) з А(р, р), (5) 
де 0- Т(р ро) - бо відносна температура тіла; аг- і коефіцієнт температу- 
Пак с» Со су 


2 
- ; К 
ропровідності Та» пк Осо Вес я (сло б 


Фа Уа 


Уа .. С(ф, 2), 8 ФІ ,Ф)є С(0,2л). 


Тоді рішення крайової задачі (2)-(5) д(р,ф. 2,1) є двічі | неперервно 


диференційованим пор і оф, 2, один раз по і в області Д і неперервним на ДР |З, 


тобто | б(рирузьг)є С(р)ПІ СР), а функції Ср, 2) Є( рр), Ар, р), руту ог) 
можуть бути розкладені в комплексний ряд Фур'є | 31: 
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(р,9. б 9, яю 
С(9, | СО | 
р | Фр) екоітрі о 
Х(р,ф) Х, (р) 
і 9, (рот) М(р.р.от) 
знане 
А (р») Х(р.ф) 


0 Мои г)-0 "ЛО РЕЯРА б Гру з, г), с ( Се вла 


9, (р)- 6 (руч ів) | Ар) Му (руч їм (р) 


З огляду на те, що др, Ф, 2, г) функція дійсна, обмежимося надалі розглядом 
9 (р, 21) для по-0,1,2,..., тому що 9 (р, 21) і 9 (рол) будуть комплексно 


спряженими |3). Підставляючи значення функцій з (6) у (2)-(5), у результаті 
одержимо систему диференціальних рівнянь: 


261) 
д 


28 здано пт, 2580. «09879. «9:80 1980 зео- 0Ю 
Ї гі? п 2 


- ЗР 
ді | В2| до? р др вит да? (7) 


з початковими умовами 


28 (р.а) 


| 9 (р, 20)-0, Ей (8) 
іграничними умовами 
а є с) (9) 
дід (ро) «вїЙ(р) | ві(рл.л)- (р) (19) 
де 30) 2 -дп; 42 зап; ту 22, т, 21; із-1,2. 
Для розв'язання крайової задачі (7)-(10) застосовуємо інтегральне перетворення: 
ПицадЗ ОЇникорьс) р Гра (п) 


Власні функції г21РРЖр р, 1) і власні значення ш,; Знаходяться із розв'язку задачі 


Штурма-Ліувілля: 


920 100 п? Я 
но дао 0 а2) 
др" РОр р 02 


О(нухр0)- 0, шу жр.1)- 0, О(р, обо а)- 0. | (3) 
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Власні функції Оших ; р, 1) і власні значення ш, Задачі Штурма-Ліувілля (12) - 


(13) знаходяться за формулами, які приведені в |4|, формула оберненого 
перетворення має вигляд: 


"р,г)- У Сини? й ГАР'НИВІ а4 
па РР 2) 


Застосовуючи до системи диференціальних рівнянь (7) інтегральне 
перетворення (11), у результаті одержуємо систему звичайних диференціальних 
рівнянь: 


2 (т,) 24 (1) - 
дт) у ав, ін а"0, -оі) -ш2,8Ф (15) 


"ДОГ РИ )- 0, БУ -0, (16) 
де 
у (5... 00(р,Є(о), 20) дОи ур, 
ом з Ї о а З) ій («Ман я До «-2--6) орі) ар 
0 р 1 02 02 


Криволінійний інтеграл обчислюється по замкненому додатно орієнтованому 


контуру (рис.2.) 
2 


р «6(2) 


0 р 
Рис.2. Замкнутий контур з твірною лінією г - 2(2) 
Застосовуємо до системи диференціальних рівнянь (15) з початковими 
умовами (16) інтегральне перетворення Лапласа (51: 


ор 


(з) Г/(т)е""" аг. 


0 


чо) 


У результаті одержуємо систему алгебраїчних рівнянь відносно 919) 
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з юн 6 
РОАРО РКО т, з0т З ПИРО ТРИ Не -819 і (7) 
И пк 
де і-1,2; 4ик ек ; 
Розв'язавши систему рівнянь (17), одержуємо: 
Зі) 60 (вуз? в Р СУ ті (інт, ) 
В - , (18) 


РРО НИ Уа? п ба ГК 
де ак о (і - 1,2.) 
к 


Застосовуючи до зображення функцій (17) формули оберненого перетворення 
Лапласа, одержуємо оригінали функцій: 


бід с ЛЬ Мо, ін, Мет вену о 
рф бе ЗВ Ману ані М бнузу і) 
і) 09) 
бука Знаю МВВ, Меоузу оуані) (руту ні 
р Уецо НВОНО, Ману рані у Мат беуну ні 


РА! 00) 


-1 
0.55 ) пк 


де а (з; )- 5 - з ЯЗначення 5) для і-1,2,3,4 визначаються за 
(дез; ч1) (г, ап) 


формулами 


(г аті а т) а з гуопі)? Ат (г аті з уж Ча чо таті) Ро 


512 - 534 - г 


ХЕ 


Я 


Таким чином, з урахуванням формул обернених перетворень (6) і (14), 
одержуємо температурне поле тіла обертання, яке обертається з постійною кутовою 
швидкістю и) навколо осі 07, з урахуванням кінцевої швидкості поширення тепла: 


ик Ра :) 


(рф. 21)- 242128 (набу і -80)(р лют) -ехр(іпф), 


і ЇоарьЇ| 
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де значення 8(ши 1) 88 1) визначаються по формулах (19),(20). 


Висновки 

У статті, за допомогою розробленого нового інтегрального перетворення, 
знайдено температурне поле ізотропного тіла, яке обертається з постійною кутовою 
швидкістю и) навколо осі 0/, з урахуванням кінцевої швидкості поширення тепла, 
у вигляді збіжних ортогональних рядів по функціях Фур'є. Знайдений розв'язок 
узагальненої крайової задачі теплообміну ізотропного тіла, яке обертається, з 
урахуванням скінченності величини швидкості поширення тепла може знайти 
застосування при модулюванні температурних полів, які виникають у багатьох 
технічних системах (супутники, сортопрокатні валки, ротори енергетичних 
агрегатів, дискові гальма та ін.). 
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КЕ5ОМЕ 

М.С. Вегдапук 

Машетабісаї тодеї апа яепегаїйтед 5о0їшкіоп теїоай Рігспеїе ргобіетя ої реаї 
ехспапге ої гадіайоп роду 

Іа Фіз рарег, Фе саїсціайоп ої Фе поп-5іайопагу іетрегаїшге Не!д ої ап і5оїгоріс 
роду ої гоіайоп уїб а Кпомп едиайоп ої Ше сгеакиге Ппе їп а суппатгіса! соогаїпаіе 
зувіет 15 сопзідегей, мрісі гокакеє мій а соп5апі апеціаг уеіосіїу агойпа Фе О7 ахіз, 
сакіпє їпіо ассоппі Ше па! гаїе ої реаг ргорабайоп. Трегтаї-ррузісаї ргорегіез ої мпіср 
до пог дерепа оп іеппрегаїиге апа іп(егпа! 50йгсе5 ої Пеає аге аб5епі. Ас Ше іпіца! 
тотепі, Ше (епірегашге ої Бе Боду 15 соп5іапі, апа оп Бе оці5іде ої (Фе суппаєг 
іетрегаїиге 15 Кпом/п апа дое5 пої дерепа оп те. 

Бог Бе Пг5: йте, а пзафетайсаї тоде! ої (ептрегаїиге дї5ітібийоп5 ої ап і5оігоріс 
роду ої гокайоп мій а Кпом/уп едиайоп ої а сгеаїиге Ппе іп а суппагіса! соогаїпаіе 
зубіет, міс гоіаїе5 мій а соп5капі апеціаг усіосіїу агоцпа Ше О7 ахі8, (аКіпє іпіо 
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ассоцпі Ше Ппа! гаїе ої реаг ргорагайоп, 15 рге5епіед а5 а Бойпдагу уаїше ргобіет ої 
таШетайсаї! рпузіс5 Їог а пурегройс ФіНегепиа! едиайоп ої іегта! сопацсйуїсу мій 
Боипдагу сопаїопз ої Ше Нг5/ 5огі ої То 50Їуе Ше Боппдаагу-уаїше ргобіет, Ше дезігед 
іетрегаїиге Пе!4 15 герге5епіеа а5 а сотпріех Боигіег 5егіе5. А пему іпіеєга! (гап5їогт Ба5 
Бееп ЧФеуеіоред Бу міс Ше (етрегаїшге Пе! ої Ше Їогт ої сопуегеепі огіпобопа! 5егіе5 
їп (егт5 ОЇ зресіа! Боигіег Гипспопея 15 Гоцпа. 

Тре 580Їийоп ої Бе єепегайоед Бойпаагу-уаїце Реаї ігап5їег ргобіет Їог ап і5оїгоріс 
госайпе Боду, мПпісП 15 гоіаїед, (акіпє їпіо ассоцпі Ше Ппігепе85 ої Ше уеіосіу ої Пеаї 
ргораєганйоп, сап Бе ц5еад Їог плодиіайоп ої гетрегаштге Пед5 Ша оссиг іп плапу (есппісаї 
зубіетя (5а хе Шіез, гоПег плій8, гоіог5 ої роу/ег абетеваїез, ді5К БгаКез, ес.) 
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